Multiwavelength Observations of the VHE Blazar 1ES 2344+514 
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ABSTRACT 

Multiwavelength observations of the high-frequency-peaked blazar 1ES 2344+514 were per- 
formed from 2007 October to 2008 January. The campaign represents the first contemporaneous 
data on the object at very high energy (VHE, E >100 GeV) 7-ray, X-ray, and UV energies. 
Observations with VERITAS in VHE 7-rays yield a strong detection of 20 with 633 excess 
events in a total exposure of 18.1 hours live-time. A strong VHE 7-ray flare on 2007 December 7 
is measured at F(>300 GeV) = (6.76 ± 0.62) x 10 -11 ph cm" 2 s _1 , corresponding to 48% of the 
Crab Nebula flux. Excluding this flaring episode, nightly variability at lower fluxes is observed 
with a time-averaged mean of F(>300 GeV) = (1.06 ± 0.09) x 10~ n ph cmT 2 s" 1 (7.6% of the 
Crab Nebula flux). The differential photon spectrum between 390 GeV and 8.3 TeV for the 
time-averaged observations excluding 2007 December 7 is well described by a power law with 
a photon index of V = 2.78 ± 0.09 sta t ± 0.15 syst . On the flaring night of 2007 December 7 the 
measured VHE 7-ray photon index was T — 2.43 ± 0.22 stat ± 0.15 syst . Over the full period of 
VERITAS observations contemporaneous X-ray and UV data were taken with Swift and RXTE. 
The measured 2-10 keV flux ranged by a factor of ^7 during the campaign. On 2007 December 8 
the highest ever observed X-ray flux from 1ES 2344+514 was measured by Swift XRT at a flux of 
F(2-10 keV) = (6.28 + 0.31) x 10~ n erg cm~ 2 s -1 . Evidence for a correlation between the X-ray 
flux and VHE 7-ray flux on nightly time-scales is indicated with a Pearson correlation coefficient 
of r = 0.60 ± 0.11. Contemporaneous spectral energy distributions (SEDs) of 1ES 2344+514 
are presented for two distinct flux states. A one-zone synchrotron self-Compton (SSC) model 
describes both SEDs using parameters consistent with previous SSC modeling of 1ES 2344+514 
from non-contemporaneous observations. 

Subject headings: galaxies: BL Lacertae objects: individual: 1ES 2344+514 
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1. Introduction 

The majority of blazars detected at very high 
energy (VHE, E >100 GeV) 7-rays are high- 
frequency-peaked BL Lac objects (HBLs), with 
currently ~25 HBLs from a total of -30 VHE 
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blazarsQ The short time-scale variability seen 
in the broadband spectral energy distributions 
(SEDs) of blazars is explained by highly rel- 
ativistic plasma jets oriented close to the line 



ativistic plas ma jets orientea close to tne line 
of sight (jBlandford fe Konigll [l979h . HBLs are 
blazars exhibiting synchrotron radiation peaking 
typically at UV to X-ray energies and a second 
SED component peaking at GeV to TeV energies. 
Leptonic models (Coppi 1992t Bottcher fc Chang] 



2002; iKrawczvnski et al.ll2004j ) describe the high- 



energy peak as inverse Compton upscattering of 
low-energy photons by electrons, while hadronic 
models attribute the emission to proton-induced 
cascades, synchrotron radiation by protons, pion 
production in dens e clumps of jet plasma, or 
curvature radiation (Mannheim 19931: Aharonian 
20001: iMiicke fc Protherodl200ll . Detailed multi- 



wavelength studies from optical to 7-ray energies 
of the temporal and spectral variability of HBLs 
promise to constrain the physical parameters of 
the underlying particle distributions, particularly 
the Doppler factor and magnetic field strength 
(|Tavecchio. Maraschi fc Ghiselliiiilll998l ). 

The HBL 1ES 2344+514 was first detected in 
the Einstein Slew Survey at X-ray ene rgies (0.2-4 
keV) and has a redshift of z = 0.044 (jElvis et al. 
Il992t IPerlman et al.lll996l) . VHE 7-ray emission 



from 1ES 2344+514 was discovered by the Whip- 
ple 10 m telescope in energies >350 GeV at a 6.0 
a detection le yel during a 1 day f lare on 1995 



December 20 ( Catanese et al. 19981 ). X-ray ob- 



servations of 1ES 2344+514 with BeppoSAX re- 
vealed variability in the 2-10 keV flux by a fac- 
tor of —2 during a 7-da y period in 1996 Decem- 
ber (|Giommi et al.ll2000h . 1ES 2344+514 was later 
observed by Swift in 2005 during April, May, and 
December at a lower flu x than in 1996 December 
( Tramacere et al 1 l2007l) . The X-ray spectrum is 
reasonably well described by an absorbed power 
law, with photon indices ranging from L w 1.8- 
2.3. Marginal evidence for optical variability in 
1ES 2344+514 is seen from the sparse data set 
of two observations in 199 8 ( Nilsson et al.l fl999; 
Falomo fc Kotilainenl 19991 ). six nights in 2000 



(|Xie et al l l2002h. and one observation on 2007 
January 12 ([Gupta et al 1l2008h . In the radio band, 
VLA observations reveal an unusual morphology 



1 The TeVCat catalog of VHE 7-rays sources is available on- 
line at: |http : / / tevcat . uchicago . edu | 
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at arcsec scales for a blazar of two compact high 
brightness emission regions connected via a dif- 
fuse halo, more reminiscent of a radio galaxy 



I Rec tor, Gabuzda fc Stock e 2003; Giroletti et al 



2004). High resolution VLB A observations show 
a well-collimated jet extending ~10 pc from the 
core, which then bends 25° and broadens into a 
cone of ~35° opening angle. In high energy (HE, 
0.1>E<300 GeV) 7-rays, 1ES 2344+514 is asso- 
ciated with the source 1FGL J2347.1+5142 from 



the F ermi LAT First Source Catalog (jAbdo et al 
2010h . detected at a significance of 10.6 a and 



showing a hard HE 7-ray spectrum of 1.57 ± 0.17 
with no indication of HE 7-ray flux variability dur- 
ing the period of 2008 August to 2009 July. 

Evidence for long-term VHE 7-ray flux variabil- 
ity in 1ES 2344+514 is seen in observations be- 
tween 1995 and 2005 by the Whipple 10 m tele- 
scope, HEGRA, and MAGIC. The Whipple 10 m 
7-ray photon spectrum on the flaring night of 1995 
December 20 is best-fit by a power law over the 
energy range 0.8-9 TeV, with a photon index Y 



2.54 ± 0.17 sta t ± 0.07 syst (|Schroedter et al.ll2005l ) 



The measured flux F(>1 TeV) = (3.3 ± 0.7) 
10~ n ) ph cm~ 2 s _1 for the flare night is a fac- 
tor of 1.7 ± 0.3 higher than the Crab Nebula flux. 
Excluding 1995 December 20, the time-averaged 
Whipple 10 m data between 1995 and 1997 yielded 
a marginal 4 a detection at a significantly lower 
flux level, c orresponding to ~10 % of the Crab 
Nebula flux (jCatanese et alJll998h . Further VHE 
7-ray observations of 1ES 2344+514 with HEGRA 
in 1997, 1998, and 2002 resulted in a 4.4 a detec- 
tion in a low flux le vel state of (3.3 + 1.0)% of the 



Crab Nebula flux (jAharonian et al.1 l2004h . Ob 



servations in 1995 with TACTIC above 1.5 TeV 
in 2004 and 2005 yielded a weakly constraining 
flux upper limi t near the Whipple 1 m detec- 



tion flux level (|Godambe et all 120071 ). In 2005 
August to 2006 January, the MAGIC telescope 
measured a time-averaged photon spectrum be- 
tween 0.14-5.4 TeV at a low flux of 5% of the 
Crab Nebula flux characterized by a power-law 
with T = 2.95 ± 0.12 stat ± 0.20 syst (I Albert et al ' 



2007). This article presents the first multiwave- 
length campaign on 1ES 2344+514, with contem- 
poraneous UV, X-ray and VHE 7-ray observations 
over a several month long period. 
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Fig. 1. — The distribution of 8 2 for 011- 
source events (points) and normalized off-source 
events (shaded region) from observations of 
1ES 2344+514. The dashed line represents the ap- 
plied cut on the 9 2 parameter. 

2. VERITAS Observations 

VERITAS is an array of four 12 m diameter 
imaging atmospheric Cherenkov telescopes lo- 
cated at the Fred Law rence Whipple Obse rvatory 



in Southern Arizona ( Weekes et al. 20021 ). Each 



VERITAS camera contains 499 pixels (0.15° di- 
ameter) , and has a field of view of 3.5° . VERITAS 
is sensitive over an energy range of ^100 GeV to 
~30 TeV with an energy resolution of 15-20%, and 
an angular resolution (68% containment ) of less 
than 0.14° per event (jAcciari et al.ll2008l ). VER- 
ITAS can detect VHE 7-ray fluxes of 5% and 1% 
of the Crab Nebula flux at the 5 a level in <2.5 
and <50 hours, respectively. 

The VERITAS observations of 1ES 2344+514 
presented here were taken on 37 nights between 
2007 October 4 and 2008 January 11. After ap- 
plying quality-selection criteria, the total exposure 
is 18.1 hours live-time. Data-quality selection re- 
quires clear atmospheric conditions, based on in- 
frared sky temperature measurements, and nor- 
mal hardware operation. During each night, the 
quality-selected data ranged from 0.3-2.7 hours 
live-time. The zenith angle of observations ranges 
from 19-48°, with a mean of 28°. All data were 
taken during moon-less periods in wobble mode 
with pointings of 0.5° from the blazar alternat- 
ing from North, South, East, and West direc- 
tions to_ejiable_S2mi^aiieous background estima- 
tion ( Aharonian et al.l l2001). 
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Data reductio n follow ed the methods described 
in Acciari et al. ( 2008|). Signals in e ach event 
are first calibrated ( Holder et al.l 2006 ). and the 
images are parameterized (jrlillasl Il985h . The 
7-ray direction and air shower impact param- 
eter on the ground are th en reconstructed us- 
ing stereoscopic techniqu es ( Hofmann et al. 1999t 
Krawczvnski et all 120061 ). The background of 
cosmic-rays is rejected with a very high efficiency 
using event-by-event cuts on the arrival direction 
(9 2 ), mean scaled width and length, integrated 
charge (size), and location of the image centroids 
in the camera (distance). The cuts applied here 
are optimized a priori for a source strength of 
10% of the Crab Nebula flux and a similar photon 
index to the Crab Nebula. The energy of each 
event is reconstructed using looku p tables from 



Mont e Carlo simulations of 7-rays (jAcciari et al 



120081 ). Results fr om two independe nt VERITAS 
analysis packages (jDaniel et al.ll2007l ) yield results 
consistent with those presented here. 

A total of 1275 on-source and 4494 off-source 
events were measured, with an on-off normaliza- 
tion of a = 0.143. This corresponds to a to- 
tal excess of 632 events from the direction of 
1ES 2344+514. The statistical significance of this 
excess is 20.2 standard deviations (a). Figure Q] 
shows the 9 2 distribution of the total data set, 
which is the squared angular distance between 
the reconstructed event direction and the nomi- 
nal source position. The shape of the excess is 
consistent with a simulated point source for VER- 
ITAS. 

3. RXTE Observations 



The X-ray satellite mission RXTE ([Bradt et al 



19931 ) observed 1ES 2344+514 between 2007 Octo- 
ber 7 and 2008 January 11 (ObsID 93132). Data 
are presented from Pr oportional Counter Unit 2 of 
the PC A instrument ( Jahoda et al.lfl996h since for 
nearly all observations the other PCUs were not in 
operation. The 52 nightly PCA observations were 
taken in snapshots ranging from 0.11-1.00 hours 
live-time, and are listed in Table [TJ Data reduc- 
tion is performed with the HEAsoft 6.5 package. 
Only the top layer (X1L and X1R) signal is used. 
The data are filtered following the standard cri- 
teria advised by the NASA Guest Observer Facil- 
ity. Background data are parameterized with the 



pcabackest tool using the most recent model for 
faint sources. The photon spectrum of each ob- 
servation is extracted using the saextrct tool. Re- 
sponse matrices are generated using pcarsp with 
the latest calibration files. 

4. Swift Observations 



Swift (|Gehrels et al.l 120041 ) Target of Oppor- 
tunity (ToO) observations of 1ES 2344+514 from 
2007 October 27 to 2008 January 1 were triggered 
by the VERITAS detection. The Swift data set 
consists of eight snapshot observations of 10-45 
minute duration each, as listed in T able [TJ All 
Swift XRT data ([Burrows et al.ll2005l ) are reduced 
using the HEAsoft 6.5 package. Event files are cal- 
ibrated and cleaned following the standard filter- 
ing criteria using the xrtpipeline task and applying 
the most recent Swift XRT calibration files. All 
data were taken in Photon Counting (PC) mode, 
with grades 0-12 selected over the energy range 
0.3-10 keV. Due to photon pile-up in the core of 
the point spread function (PSF) at count rates 
>0.5 counts s _1 (PC mode), the source events are 
extracted from an annular region with inner ra- 
dius ranging from 2-6 pixels and an outer radius 
of 30 pixels (47.2 arcsec). Background counts are 
extracted from a 40 pixel radius circle in a source- 
free region. Ancillary response files are generated 
using the xrtmkarf task, with corrections applied 
for the PSF losses and CCD defects. The latest 
response matrix from the XRT calibration files is 
applied. To ensure valid \ 2 minimization statistics 
during spectral fitting, the extracted XRT energy 
spectra are rebinned to contain a minimum of 30 
counts in each bin. 

UVOT observations were taken in the pho- 
tometric bands o f UVW1, UVM2, and UVW2 
( Poole et al. 2008). The uvotsource tool is used 



to extract counts, correct for coincidence losses, 
apply background subtraction, and calculate the 
source flux. The standard 5 arcsec radius source 
aperture is used, with a 20 arcsec background re- 
gion. The source fluxes are de reddened us i ng th e 
interstellar extinction curve in iFitzpatric k (1999). 
In the UV filters only a low level of host galaxy flux 
is evident ([Tramacere et al.l 120071 ). so no correc- 
tions are needed. The uncertainty in point spread 
function variations with time is fixed for this anal- 
ysis at 15%. Observations were taken with just one 
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Fig. 2. — Differential photon spectrum of VHE 7- 
rays for 1ES 2344+514. The circles represent the 
time-averaged VERITAS data from 2007 October 
to 2008 January, excluding the bright flare on 2007 
December 7. The open squares represent VERI- 
TAS data from the flare night. Shown for compar- 
ison are the 68% confidence intervals for the best- 
fit power law models to 1ES 2344+514 data from 
MAGIC (dark grey band) in 2005 August to 2006 
January (j Albert et al.ll2007l ). and from the Whip- 
ple 10 m (light grey ba nd) during a high flux s tate 
on 1995 December 20 llSchroedter et al.ll2005l ). 
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Fig. 3.— X-ray spectrum of 1ES 2344+514 from 
Swift XRT and RXTE PCA observations. Shown 
in the top panel are example moderate-flux and 
high-flux spectra identified in the plot by instru- 
ment and observation date. The differential rates 
are scaled in three cases for viewing purposes. For 
each energy bin, the statistical significance is >3 a, 
and the horizontal line represents the best-fit for- 
ward folded absorbed power law model. The bot- 
tom panel shows the ratio of the data and model 
values. 



UV filter during each pointing, and only UVM2 
data were taken on multiple nights. Within the 
conservative errors, no variability is seen between 
the four UVM2 filter observations. 

5. Spectral Analysis 

The VERITAS differential photon spectrum 
over the energy range of 390 GeV to 8.3 TeV for 
the time-averaged data excluding a strong flaring 
night on 2007 December 7 is shown in Figure [2] A 
power law model fit to the data is best described 
by the form dN/dE = F • (E/lTcV)~ r yielding 
a % 2 /dof of 8.39/4 for a flux normalization con- 
stant of F = (3.38 ± 0.23 sta t ± 0.70 syst ) x 10~ 12 
ph cm~ 2 s _1 TeV" 1 and a photon index of T = 
2.78±0.09 s tat±0.15 sys t. A fit to a power law with 
exponential cutoff model, which contains one ad- 
ditional degree of freedom, does not yield a signifi- 
cantly better fit. On 2007 December 7 the weather 
conditions were affected by moving cloud coverage, 
so the observations do not pass the standard data- 
quality criteria. However, online analysis showed 



that the source was in the most active state yet 
observed. We therefore include these data and 
increase the systematic uncertainties in flux and 
spectral shape. Conservative systematic uncer- 
tainties in the flux of 30% and photon index of 0.2 
are determined, based on the results of observa- 
tions of the Crab Nebula taken in similarly vari- 
able conditions. The VHE 7-ray spectrum from 
the high flux night is well described by a power law 
(x 2 /dof of 1.27/3) with a slightly harder photon 
index T — 2.43±0.22 sta t±0.2 sys t, and a higher flux 
normalization ofF = (17.3±1.92 s t a t±5.19 sys t) X 
10" 12 ph cm" 2 s- 1 TeV- 1 . 

The best-fit results presented here are tested 
with a dedicated Monte-Carlo study, which ex- 
plores the uncertainty ranges for the two measured 
spectra. In the Monte-Carlo study photons are se- 
lected from a power-law distribution matching the 
measured distribution convolved with the effective 
area as a function of energy. The selected photon 
energies are then smeared with a Gaussian cor- 
responding to the energy resolution of ~15% at 
these energies. The number of photons selected is 
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Tabic 1: Best-fit spectral parameters for the Swift XRT and RXTE PC A data. 



Date 


Start 


Exp. 


r 


K-lkcV 


Xr/doi 


F(2-10 kcV) 




(hr:min UT) 


(kscc) 




(10~ 2 ph cm" 2 s" 1 keV -1 ) 




(10 — 11 erg cm -2 s — 1 ) 


Swift XRT 














2007-10-27 


10:24 


2.66 


2.33 ± 0.05 


0.62 ± 0.03 


1.11/36 


0.96 ± 0.06 


2007-11-03 


03:02 


1.97 


2.32 ± 0.06 


0.87 ± 0.04 


1.02/27 


1.37 ± 0.09 


2007-11-10 


05:03 


0.63 


2.26 ± 0.13 


1.26 ± 0.10 


0.69/8 


2.18 ± 0.31 


2007-11-17 


02:28 


1.65 


2.16 ± 0.06 


0.76 ± 0.04 


0.71/23 


1.51 ± 0.14 


2007-11-30 


03:46 


1.85 


2.10 ± 0.07 


0.91 ± 0.05 


1.51/23 


1.99 ± 0.14 


2007-12-08 


03:20 


2.22 


1.97 ± 0.05 


2.37 ± 0.10 


0.77/33 


6.28 ± 0.31 


2007-12-22 


04:28 


2.23 


2.18 ± 0.05 


1.00 ± 0.04 


1.12/36 


1.93 ± 0.13 


2008-01-01 


00:35 


1.37 


2.10 ± 0.08 


1.26 ± 0.07 


1.07/18 


2.72 ± 0.17 


RXTE PCA 














2007-10-07 


06:03 


2.56 


2.48 ± 0.26 


0.75 ± 0.35 


0.52/38 


0.95 ± 0.26 


2007-10-08 


06:56 


2.88 


2.65 ± 0.20 


1.24 ± 0.38 


0.58/38 


1.24 ± 0.15 


2007-10-09 


06:14 


2.80 


2.31 ± 0.13 


1.04 ± 0.23 


0.54/38 


1.66 ± 0.10 


2007-10-10 


05:53 


2.67 


2.17 ± 0.12 


0.89 ± 0.19 


0.47/38 


1.76 ± 0.14 


2007-10-11 


03:42 


2.16 


2.57 ± 0.15 


1.65 ± 0.40 


0.42/38 


1.85 ± 0.11 


2007-10-12 


03:14 


0.21 


2.57 ± 0.19 


1.34 ± 0.41 


0.66/38 


1.48 ± 0.15 


2007-10-13 


04:52 


2.02 


2.26 ± 0.13 


1.11 ± 0.24 


0.63/38 


1.92 ± 0.15 


2007-10-14 


02:17 


1.86 


2.39 ± 0.12 


1.63 ± 0.33 


0.62/38 


2.35 ± 0.14 


2007-10-15 


03:56 


1.50 


2.49 ± 0.16 


1.61 ± 0.41 


0.53/38 


2.02 ± 0.20 


2007-10-16 


04:35 


1.84 


2.25 ± 0.18 


0.80 ± 0.27 


0.53/38 


1.41 ± 0.24 


2007-10-19 


07:09 


3.62 


2.42 ± 0.14 


1.02 ± 0.23 


0.51/38 


1.39 ± 0.13 


2007-10-20 


08:18 


3.15 


2.61 ± 0.15 


1.51 ± 0.35 


0.59/38 


1.59 ± 0.12 


2007-10-21 


07:57 


3.49 


2.36 ± 0.16 


0.81 ± 0.23 


0.56/38 


1.22 ± 0.11 


2007-11-01 


03:06 


2.59 


2.19 ± 0.15 


0.73 ± 0.21 


0.79/38 


1.40 ± 0.13 


2007-11-02 


05:53 


1.39 


2.41 ± 0.18 


1.25 ± 0.37 


0.62/38 


1.73 ± 0.16 


2007-11-03 


07:06 


2.37 


2.72 ± 0.19 


1.65 ± 0.48 


0.53/38 


1.50 ± 0.14 


2007-11-04 


04:56 


2.85 


2.19 ± 0.14 


0.75 ± 0.20 


0.50/38 


1.43 ± 0.11 


2007-11-05 


06:35 


0.80 


2.20 ± 0.29 


0.68 ± 0.40 


0.64/38 


1.28 ± 0.45 


2007-11-06 


04:02 


2.48 


2.23 ± 0.11 


1.10 ± 0.21 


0.82/38 


2.00 ± 0.11 


2007-11-07 


06:36 


0.31 


2.43 ± 0.12 


1.28 ± 0.25 


0.49/38 


1.75 ± 0.09 


2007-11-08 


07:52 


2.96 


2.20 ± 0.10 


1.10 ± 0.18 


0.57/38 


2.07 ± 0.11 


2007-11-09 


07:35 


2.40 


2.10 ± 0.09 


1.16 ± 0.17 


0.59/38 


2.53 ± 0.12 


2007-11-10 


03:50 


2.53 


2.13 ± 0.08 


1.41 ± 0.18 


0.33/38 


2.94 ± 0.09 


2007-11-11 


01:49 


2.27 


2.16 ± 0.08 


1.50 ± 0.20 


0.63/38 


3.02 ± 0.13 


2007-11-12 


04:31 


2.51 


2.23 ± 0.13 


0.91 ± 0.21 


0.54/38 


1.65 ± 0.16 


2007-11-13 


02:30 


2.30 


2.11 ± 0.09 


1.19 ± 0.18 


0.71/38 


2.57 ± 0.13 


2007-11-14 


03:37 


2.40 


2.22 ± 0.11 


1.15 ± 0.21 


0.55/38 


2.12 ± 0.10 


2007-11-15 


03:10 


2.34 


2.50 ± 0.16 


1.31 ± 0.33 


0.60/38 


1.61 ± 0.15 


2007-11-30 


02:20 


0.37 


2.44 ± 0.26 


1.91 ± 0.81 


0.74/38 


2.54 ± 0.37 


2007-12-01 


01:51 


0.50 


2.59 ± 0.33 


1.62 ± 0.84 


0.38/38 


1.76 ± 0.46 


2007-12-02 


05:03 


0.58 


2.23 ± 0.18 


1.48 ± 0.43 


0.72/38 


2.65 ± 0.25 


2007-12-04 


02:31 


0.40 


2.10 ± 0.19 


1.32 ± 0.42 


0.61/38 


2.87 ± 0.39 


2007-12-05 


02:11 


1.79 


2.01 ± 0.05 


2.15 ± 0.19 


0.81/38 


5.37 ± 0.15 


2007-12-06 


03:28 


1.57 


2.06 ± 0.06 


1.99 ± 0.22 


0.69/38 


4.63 ± 0.11 


2007-12-10 


02:18 


2.03 


1.94 ± 0.06 


1.39 ± 0.15 


0.90/38 


3.90 ± 0.10 


2007-12-11 


02:50 


1.82 


1.87 ± 0.04 


1.97 ± 0.15 


0.62/38 


6.10 ± 0.11 


2007-12-13 


03:38 


3.39 


1.86 ± 0.04 


1.54 ± 0.10 


0.67/38 


4.85 ± 0.07 


2007-12-12 


05:41 


1.23 


1.97 ± 0.06 


1.96 ± 0.21 


0.46/38 


5.20 ± 0.13 


2007-12-14 


04:35 


3.33 


2.08 ± 0.06 


1.53 ± 0.15 


0.65/38 


3.42 ± 0.09 


2007-12-15 


04:53 


0.86 


2.17 ± 0.13 


1.57 ± 0.33 


0.91/38 


3.10 ± 0.23 


2007-12-28 


03:46 


1.06 


2.14 ± 0.11 


1.48 ± 0.28 


0.72/38 


3.06 ± 0.17 


2007-12-29 


04:53 


1.17 


1.96 ± 0.15 


0.69 ± 0.22 


0.47/38 


1.86 ± 0.23 


2007-12-30 


04:25 


0.88 


1.91 ± 0.09 


1.40 ± 0.22 


0.56/38 


4.04 ± 0.22 


2008-01-01 


01:52 


0.86 


2.06 ± 0.11 


1.43 ± 0.27 


0.85/38 


3.30 ± 0.15 


2008-01-02 


02:26 


2.66 


2.12 ± 0.07 


1.40 ± 0.17 


0.53/38 


2.96 ± 0.11 


2008-01-03 


03:34 


2.91 


1.97 ± 0.06 


1.23 ± 0.13 


0.67/38 


3.26 ± 0.10 


2008-01-04 


01:53 


1.71 


2.06 ± 0.13 


0.87 ± 0.20 


0.65/38 


2.03 ± 0.12 


2008-01-05 


02:35 


2.83 


1.97 ± 0.07 


1.02 ± 0.13 


0.51/38 


2.71 ± 0.11 


2008-01-06 


02:08 


2.75 


2.31 ± 0.12 


1.08 ± 0.22 


0.41/38 


1.76 ± 0.11 


2008-01-07 


03:48 


1.98 


2.33 ± 0.14 


1.16 ± 0.27 


0.68/38 


1.83 ± 0.12 


2008-01-10 


02:23 


1.50 


2.25 ± 0.10 


1.78 ± 0.29 


0.52/38 


3.12 ± 0.14 


2008-01-11 


03:24 


2.77 


2.12 ± 0.07 


1.53 ± 0.17 


0.74/38 


3.25 ± 0.10 



matched to the observed number of photons with 
a Poisson scatter. Five million synthetic photon 
spectra are created using similar binning to the 
original spectra, and the best-fit power law index 
and normalization parameters are extracted for 
each synthetic spectrum providing dense coverage 
of the parameter space. The Monte-Carlo mean 
value and one dimensional standard deviation of 
the power law parameters are consistent with the 
best-fit parameters of the measured spectra. For 
the time-averaged spectrum excluding the flaring 
night of 2007 December 7, the Monte-Carlo study 
yields a mean value and 1-er error for the flux nor- 
malization of F ,mc = (3.35 ± 0.30) x 10~ 12 ph 
cm~ 2 s _1 TcV -1 and for the the photon index of 
T MC = 2.76 ± 0.15. For the high flux 2007 De- 
cember 7 spectrum, the Monte-Carlo study yields 
F ,mc = (17.3±2.67) x 10~ 12 ph cm" 2 s" 1 TeV" 1 
and r M c = 2.40 ± 0.28. 

Previous VHE 7-ray observations of 1ES 2344+514 
with the MAGIC telescope between 2005 August 
and 2006 January measured a constant flux at a 
slightly lower level than the time-averaged (flare 
excluded) flux measured here, and a consistent 
photon index of T = 2.95 ± 0.12 sta t. ± 0.2 syst . 



pared to a joint fit with a fixed G alactic column 



densi ty N H ,Gai of 1.50 ■ 10 21 cm" 2 (jKalberla et al 



(| Albert et all 120071 ). The high flux VHE 7-ray 
photon spectrum measured by the Whipple 10 
m on 1995 December 20 with a photon index of 
r = 2.54 ± 0.17 sta t ± 0.07 syst (jSchroedter et al ' 



2005) is marginally harder than the low flux spec- 
trum and agrees with the high flux spectrum mea- 
sured here. The VHE observed 7-ray spectral 
shape and intensity are modified by absorption 
from pair production interactions with the in- 
frared compo nent of the extragalactic background 
light (EBL) (|Gould fc Schrederlll967l) . Using the 
redshift z = 0.044 for 1ES 2344+514, we calculate 
a d e-absorbed spectrum base d on the EBL model 
of iFranceschini et aL ( 2008 ). yielding intrinsic 
photon indices for the VERITAS low and high 
flux states of r int j ow « 2.5 and r inti high ~ 2.1, 
respectively. 

X-ray spectral analysis of the Swift XRT and 
RXTE PCA data is performed with XSPEC 12.4. 
An absorbed power law model, including the phabs 
model for the photoelectric absorption, is fit to 
each spectrum. First, a joint fit of the eight Swift 
XRT over the energy range 0.4-10 keV using a tied 
column density Nh and each spectrum having a 
varying photon index and normalization is com- 



120051 ). The joint fit with tied N H yields a best- fit 
Njj of (2.06+0. 13)-10 21 cm" 2 , with a reduced X 2 of 
1.03 for 203 degrees of freedom (dof) and chance 
probability for larger chi-squared statistics (null 
hypothesis) of 37%. The alternative joint fit with 
the fixed Galactic column density Nn,Gai yields a 
higher reduced x 2 of 1.13 for 204 dof and null hy- 
pothesis of 9.2%. The best-fit Nh is then used 
in the absorbed power law model for the RXTE 
PCA spectra, which are mostly unaffected by ab- 
sorption due to the 3-20 keV energy range, and 
for all Swift XRT results presented here. Figure |3] 
shows Swift XRT (0.4-10 keV) and RXTE PCA 
(3-20 keV) spectra from example moderate-flux 
and high-flux observations. The ratio in each en- 
ergy bin of the data to the absorbed power law 
model is shown for the example spectra, with the 
full list of reduced % 2 values per dof for each obser- 
vation listed in Table [1] The model describes well 
the spectra, yielding an average chance probability 
for larger chi-squared statistics of 37% and 51% for 
the Swift XRT and RXTE PCA data, respectively. 
The measured Swift XRT 0.4-10 keV photon in- 
dices range from T = 2.33 + 0.05 to T = 1.97 + 0.05 
for observations differing in 2-10 keV integrated 
flux by a factor of 6.54 ± 0.52. These results are 
consistent with the BeppoSAX 0.5-10 keV pho- 
ton indices of T w 1.8-2.3 at similar flux levels for 
an absorbed po wer law model with fixed Galactic 
column density ( Giommi et al.l l2000). Detailed re- 
sults from the nightly RXTE PCA and Swift XRT 
spectra focusing on the connection between the 
significant X-ray spectral and flux variability are 
presented in the following section. 

6. Light Curves 

The nightly VHE 7-ray and X-ray light curve 
of 1ES 2344+514 from VERITAS, RXTE PCA, 
and Swift XRT is shown in Figure |4] A strong 
VHE 7-ray flare is seen on 2007 December 7 
(54441.12 MJD) at an integrated flux F(>300 
GeV) = (6.8 + 0.6) x 10~ n cm~ 2 s -1 , correspond- 
ing to 48% of the Crab Nebula flux. The measured 
increase in flux of a factor of 1.9 ± 0.5 between 
the previous night and the flare night shows the 
first clear evidence of ^day time-scale VHE 7-ray 
variability from 1ES 2344+514 since the Whipple 
10 m discovery of VHE 7-ray emission in 1995 



7 



Oct. 2007 



Nov. 2007 



+ ++++ 



Dec. 2007 



Jan. 2008 
VERITAS 



370 375 380 385 390 395 400 405 410 415 420 425 430 435 440 445 450 455 460 465 470 475 480 

Date [ MJD - 54000 ] 



• RXTE PCA 
□ Swift XRT 



♦I* 

1 



♦ : • 



370 375 380 385 390 395 400 405 410 415 420 425 430 435 440 445 450 455 460 465 470 475 480 

Date [ MJD - 54000 ] 




370 375 380 385 390 395 400 405 410 415 420 425 430 435 440 445 450 455 460 465 470 475 480 

Date [ MJD - 54000 ] 



Fig. 4.— VHE 7-ray and X-ray light curve of 1ES 2344+514. Shown in the top panel are VERITAS F(>300 
GeV) nightly fluxes. The middle panel shows the 2-10 keV fluxes from observations with RXTE PCA 
(circles) and Swift XRT (open squares). The bottom panel shows X-ray photon indices from the best-fit 
absorbed power law model. 
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Fig. 5. — X-ray photon index versus 2-10 keV 
flux from RXTE PCA (circles) and Swift XRT 
(open squares) data on a log-log scale. The best- 
fit power law for the RXTE PCA (3-20 keV) data 
is represented by a solid line, and for Swift XRT 
(0.4-10 keV) with a dashed line. 
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Fig. 6.— VERITAS 7-ray flux F(>300 GeV) ver- 
sus X-ray 2-10 keV flux from nights with RXTE 
PCA (circles) and Swift XRT (open squares) data. 
A linear fit to the data is represented by the solid 
line. 



(jCatanese et al.lll998l ). Excluding the 2007 De- 
cember 7 flaring event, the average F(E>300 GeV) 
is (1.1 ±0.1) x 10 -11 ph cm~ 2 s _1 , corresponding 
to 7.6% of the Crab Nebula flux. A fit to a con- 
stant flux excluding the flaring night is rejected 
with a chance probability of 7.2 x 10~ 8 , indicating 
significant low-level variability. A measure of the 
integrated level of flux variability is the fractional 
root-mean-square (rm s) variability amplitude F var 
(|Vaughan et al.ll2003h . For the full VERITAS data 
set of nightly exposures a high level of variability 
F V ar = (75 ± 10)% is implied. Excluding the flare 
night, a still significant F var = (34 ± 16)% is de- 
termined. Searches for short-term flux variations 
within each of the nightly observations showed no 
significant variability. 

Figure 2] (lower panels) shows the 2-10 keV 
flux and photon index T measured over 3-20 keV 
from RXTE PCA and 0.4-10 keV from Swift XRT 
data. The X-ray flux is shown to be highly variable 
throughout the campaign, with F var = (51 ± 1)%. 
In 2007 December, large amplitude flaring is ev- 
ident with flux doubling time-scales of ~1 day. 
A 2-10 keV flux of (6.3 ± 0.3) x 10" 11 erg cm~ 2 
s" 1 is seen from the Swift XRT data on 2007 De- 
cember 8, representing the highest X-ray flux ever 
measured for 1ES 2344+514. Analysis of all sub- 
sequent X-ray data of 1ES 2344+514, which cur- 
rently consists solely of Swift XRT observations, 
show a 2-10 keV flux consistent with the lowest 
flux measurements presented here. Detailed re- 
sults from the more recent Swift XRT data set 
are beyond the scope of this publication, but au- 
tomated count rate light curves are publicly avail- 
able| 

X-ray spectral variability between flaring nights 
is investigated by plotting the X-ray photon in- 
dex r versus the 2-10 keV flux, shown in Figure 
[5J A logarithmic correlation of decreasing pho- 
ton index with increasing flux is suggested from a 
power law fit to the RXTE PCA and Swift XRT 
data with a x 2 /dof probability of 0.14 and 0.70 
compared to 7.9 xl0~ 20 and 3.3 xl0~ 4 for a con- 
stant photon index, respectively. Due to the lack 
of significant curvature in the X-ray spectrum, the 
peak energy E p for a vF(v) SED representation is 



2 The Swift XRT Monitoring of Fermi LAT 
Sources of Interest is available online at: 
http://www.swift.psu.edu/monitoring/ 



9 



largely unconstrained. A comparison of the pho- 
ton indices from the Swift XRT and RXTE PCA 
spectra restricted to the overlapping energy range 
of 3-10 keV is limited by the low number of en- 
ergy bins above 3 keV, so a study of the systemat- 
ically higher RXTE PCA (3-20 keV) photon in- 
dices compared with the Swift XRT (0.4-10 keV) 
photon indices at similar flux levels is inconclu- 
sive. Furthermore, the lack of purely simultaneous 
RXTE PCA and Swift XRT data does not allow 
for a detailed study of joint fits from 0.4-20 keV. 
Even during the bright flaring events in 2007 De- 
cember, the RXTE PCA spectra from 3-20 keV 
show no sign of curvature, with measured photon 
indices of ~1.9. For these high-flux states the im- 
plied peak energy E p < 10 keV agrees with the 
estimated peak energ ies derived from Be ppoSAX 
observations in 1998 ( Giommi et al. I l2000h . 

Figure EO shows the VERITAS flux F(>300 
GeV) versus RXTE PCA and Swift XRT 2-10 keV 
fluxes for nights with observations in both energy 
bands. A Pearson coefficient of r = 0.60 ±0.11 is 
calculated for the VHE 7-ray to X-ray flux points, 
suggestive of correlated variability. The best-fit 
linear model yields a slope of 0.56 ± 0.08 with a 
fit probability of 0.14. For the highest measured 
X-ray fluxes in 2007 December the weather condi- 
tions at the VERITAS site were poor, which ex- 
cludes these nightly fluxes from figure |6] Due to 
the sparse data sampling, a search for lags between 
the VHE 7-ray and X-ray emission is not investi- 
gated here. 

7. Discussion 

The broadband spectral energy distribution 
(SED) of high-frequency-peaked BL Lac objects 
(HBLs) is often modeled within a synchrotron 
self-Compton (SSC) framework. Simple one-zone 
SSC models typically predict correlated X-ray and 
VHE 7-ray variability from the underlying elec- 
tron distribution. Figure [7] shows a one-zon e SSC 
model calculation (iKrawczvnski et al. 2004) to the 
SED of 1ES 2344+514 for illustration purposes in 
two flux states. A high flux SED is constructed 
from the VERITAS photon spectrum in 2007 De- 
cember 7 and Swift UVOT and XRT data on 2007 
December 8. The time- averaged (flare excluded) 
VERITAS data is combined with Swift and RXTE 
PCA spectra on 2008 January 1 when the object 



was at a moderate X-ray flux state with respect 
to the full campaign. Shown also in Figure [7] arc 
the non-contemporaneous MAGIC data from 2005 
in a low flux state and the flaring Whipple 10 m 
data from 1995 December 20. Also shown for 
reference is the Fermi LAT spectrum from 2008 
August to 2009 July in a non-variable low flux 
state. Each model curve for the broadband data 
contemporaneous with the VERITAS data repre- 
sents the intrinsic source radiation, and all VHE 
7-ray spectra are corrected for absorption by the 



extragalactic background light (IFranceschini et al 
2008). The SSC model input parameters for the 



low (high, when different) states are: a Doppler 
factor 8 = 13 (20), magnetic field = 0.09 G (0.03 
G), emission radius = 10 16 cm, electron density = 
0.05 ergs cm~ 3 , and a broken power law electron 
spectrum with low-energies between 10 s eV and 
lQil.3 eV (10 11 - 4 eV) with spectral index m = 2.5 
(2.3), and above the cutoff energy with index ni = 
3.2 and highest energy of 10 12 eV. 

The light crossing time r = R/(c x 8) as 5 
hours for the high flux state with a Doppler factor 
8 = 20 is consistent with the ~hour time-scale X- 
ray variability seen from 1E S 2344+514 with Bep - 
poSAX on 1996 December 7 (|Giommi et alj|200d) . 
Further hourly X-ray variability in 1ES 2344+514 
could not be confirmed in this work due to the 
short (<hour) exposures with Swift and RXTE 
and moderate event statistics. In addition to the 
low and high flux state modeling from data in 
2007 October to 2008 January, the same SSC 
model was applied to non-contemporaneous low 
and high flux states derived from BeppoSAX. 



MAG IC, and Whipple 10 m spectra ([Albert et al 
2007). A reasonable agreement between model 
parameters, such as Doppler factors ~10-20 and 
a magnetic field ^0.03-0.3 G, is evident for both 
the low and high flux states with the archival 
1ES 2344+514 data, and compared with one-zone 
SSC modeling of similar moderat e-flux level flares 
in the HBLs 1ES 1959+650 dTagliaferri et al ' 



20081) . 1ES 08 06+524 (lAcciari et all l2009h. and 



1ES 1101-213 (jAharonian et al J 120071 ). A detailed 
study of the SSC model parameter sp ace (e.g. 



(|Tavecchio. Maraschi fc Ghisellinil 119981 )) is not 
pursued here due to sparseness in the SED cover- 
age and degeneracy between model parameters. 
The simple one-zone SSC modeling presented 
here adequately describes the data, however fur- 
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Fig. 7. — Spectral energy distribution of 1ES 2344+514. The high flux spectrum (triangles) are from VER- 
ITAS data on 2007 December 7 and Swift UVOT and XRT data on 2007 December 8. The time- averaged 
VERITAS spectrum and moderate flux level X-ray spectrum from Swift XRT and RXTE PCA data on 2008 
January 1 are represented by circles. An example low flux X-ray spectrum (open circles) is from 2007 Novem- 
ber 3. Non-contemporaneous VHE 7-ray spectra are shown from MAGIC data in 2005 August-December 
(open squares) and from Whipple 10 m data on the flare night of 1995 December 20 (o pen triangles). A ll 



VHE 7-ray spectra are corrected for absorption by the extragalactic bac kground light ( Franccs chini et al 
2008). Archival millimeter fluxes (open diamonds) are from SCUBA data ( Stevens fc Gear 1999). The non- 



conte mporaneous Fermi LAT spectrum from the 1FGL catalog is represented by open circles (jAbdo et al 
2010). The broadband curves are from synchrotron self-Compton (SSC) modeling to the contemporaneous 
data (jKrawczvnski et al1l2004 ). 
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ther optical, X-ray, and VHE 7-ray observations 
of 1ES 2344+514 simultaneous with Fermi LAT 
promise to offer stronger constraints to emission 
models. 

8. Conclusion 



In this paper, a ~4 month long VHE 7-ray and 
X-ray observing campaign of 1ES 2344+514 re- 
vealed significant flux variability on ~1 day time- 
scales. In particular, flux doubling between nights 
is evident in both the X-ray and VHE 7-ray bands 
in 2007 December, with the highest ever X-ray 
flux measured from 1ES 2344+514 on 2007 De- 
cember 8. Evidence of correlated nightly X-ray 
flux and VHE 7-ray flux variability is shown, as 
;enerally found in the we ll-studied HBLs Mrk 501 
Krawczvnski et al.|[2000h . Mr k 421 dFossati et al 
2008), and 1ES 1959+650 (|Krawczvnski et al 



2004), however so-called orphan flares of increased 
flux in one energy band and not the other have 
also been observed (e.g. ( Btazeiowski et al. 20051: 



Krawczvnski et al.l 120041) ) . For these campaigns 



sparse temporal sampling and insufficient VHE 
sensitivity to short time-scale variability limited 
any statements on a linear or quadratic relation- 
ship between VHE and X-ray flux expected under 
various SSC scenarios dur ing rising or decaying 
flares (|Fossati et al.l 120081 ). The 1ES 2344+514 
data presented here represents the first contempo- 
raneous multiwavelength campaign of the object 
from UV to VHE 7-rays, but is similarly restricted 
by sparse temporal sampling and broadband cov- 
erage in the spectral energy distributions. Ap- 
plying a simple one-zone SSC model reasonably 
describes the data, and is consistent with past 
modeling of the object and other HBLs. Further 
optical, X-ray, and VHE 7-ray observations of 
1ES 2344+514 simultaneous with Fermi LAT in 
the HE 7-ray band promise to offer stronger con- 
straints to emission models and high sensitivity 
broadband measurements of high-flux state flar- 
ing events. 

This research is supported by grants from the 
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